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ABSTRACT 

It is currently unknown whether the Universe was reionized by quasars or stars at z > 5. Wc point out 
that quasars can be best distinguished from steUar systems by their X-ray emission. Based on a simple 
hierarchical CDM model, we predict the number counts and X-ray fluxes of quasars at high redshifts. 
The model is consistent with available data on the luminosity function of high-redshift quasars in the 
optical and soft X-ray bands. The cumulative contribution of faint, undetected quasars in our model 
is consistent with the unresolved fraction of the X-ray background. We find that the Chandra X-ray 
Observatory might detect ~ 10^ quasars from redshifts z ^ 5 per its 17' x 17' field of view at the flux 
threshold of ~ 2 x 10~^^ erg s~^ cm~^. The redshifts of these faint point-sources could be identified by 
follow-up infrared observations from the ground or with the Next Generation Space Telescope. 

Subject headings: cosmology: theory - quasars: general - black hole physics 



1. INTRODUCTION 

Cold dark matter (CDM) cosmologies predict that the 
first baryonic objects appeared near the Jeans mass ('^ 
10^ M0) at redshifts as high as z ~ 30, and larger objects 
assembled later (Haiman & Loeb 1999b, and references 
therein). It is natural to identify these objects as the sites 
where the first stars or quasar black holes formed. Ob- 
servationally, bright quasars are currently detected out to 
z ~ 5 (Fan et al. 1999), and galaxies out to z ~ 6.7 
(Chen et al. 1999; Weymann et al. 1998). Although the 
abundance of optically and radio bright quasars declines 
at z ^ 2.5 (Schmidt et al. 1995; Pei 1995; Shaver et al. 
1996), simple models based on the Press-Schechter for- 
malism reproduce this decline and simultaneously predict 
a population of low-luminosity quasars that are too faint 
to be detected in current surveys (Haiman & Loeb 1998; 
Haehnelt et al. 1998). Preliminary evidence for such a 
population might already be indicated by the X-ray lu- 
minosity function (LF) of quasars recently measured by 
ROSAT (Miyaji et al. 1998a). The X-ray data probes 
fainter quasars than the optical data does, and has not re- 
vealed a decline in the abundance of high-redshift quasars 
as found in the optical (e.g. Schmidt et al. 1995). 

The current census of high-redshift galaxies provides an 
estimate of the global evolution of the star formation rate 
(SFR) in the redshift range < z ^ 5 (Madau 1999). 
Nevertheless, there is still considerable uncertainty about 
the evolution of the SFR at the highest redshifts probed 
(z ^ 3). Recent sub-mm observations from SCUBA in- 
dicate the existence of a substantial population of dust- 
obscured star-forming galaxies that could raise the implied 
SFR at high redshift (Barger et al. 1999, and references 
therein). This uncertainty is particularly important, since 
it leaves open the question of whether the observed popula- 
tion of stars produce sufficient ionizing radiation by z 5 



to reionize the intergalactic medium (IGM). The indica- 
tion from the present data is that the SFR either declines 
or remains roughly constant at z <; 3; however, in order 
to satisfy the Gunn-Peterson constraint for the ionization 
of the IGM, the SFR needs to rise at high redshifts above 
its value at z ~ 3 (Madau 1999). 

The presence of a substantial population of faint quasars 
at z 5 is therefore consistent with current data, and is 
suggested by cosmological models for hierarchical struc- 
ture formation. If future data would indicate a non- 
increasing cosmic SFR at z > 3, then quasars would be the 
likely sources of reionization. Recent theoretical models 
have focused on the properties of the quasar population in 
the optical (Haehnelt & Rees 1993; Haiman & Loeb 1998; 
Haehnelt et al. 1998), and in the infrared (e.g. Sanders 
1999; Almaini et al. 1999). These models were shown to 
be only mildly constrained by the lack of z ^ 3.5 quasars 
in the Hubble Deep Field (HDF; see Haiman, Madau & 
Loeb 1999). 

In this Letter we point out that quasars can be best 
distinguished from star forming galaxies at high redshifts 
by their X-ray emission. First, we illustrate the agree- 
ment between our simplest hierarchical model for quasars 
and existing ROSAT data on the X-ray LF (Miyaji et al. 
1998a) and the soft X-ray background (XRB, Miyaji et 
al. 1998b). We then show that forthcoming X-ray observa- 
tions with the Chandra X-ray Observatory {CXO; formerly 
known as AXAF) will be able to probe the abundance of 
quasars during the reionization of the Universe at z ^ 5. 



2. MODEL DESCRIPTION AND QUASAR SPECTRA 

Our model is based on the Press-Schechter mass func- 
tion of CDM halos (see Haiman & Loeb 1998 for complete 
details). We assume that each halo of mass M^aio forms a 
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central black hole of mass Afbh- The black hole mass frac- 
tion r = Mbh/Mhaio is assumed to obey a log-Gaussian 
probability distribution 



p{r) = exp[-(logr - logro)^/2cr^ 



(1) 



with logro — —3.5 and a — 0.5 (Haiman & Loeb 1999b) 
These values roughly reflect the distribution of black hole 
to bulge mass ratios found in a sample of 36 local galaxies 
(Magorrian et al. 1998) for a baryonic mass fraction of 
^ (ilb/fio) ~ 0.1. We further postulate that each black 
hole emits a time-dependent bolometric luminosity in pro- 
portion to its mass, Ln= Mi,hfq — AfbhiEdd exp(— t/to), 
where LEdd = 1-5 x 10^* Mbh/Mg erg s"^ is the Edding- 
ton luminosity, t is the time elapsed since the formation of 
the black hole, and to = lO^yr is the characteristic quasar 
lifetime. Finally, we assume that the shape of the emitted 
spectrum follows the mean spectrum of the quasar sample 
in Elvis et al. (1994) up to a photon energy of 10 keV. 
We extrapolate the spectrum up to ~ 50 keV, assuming a 
spectral slope of a=0 (or a photon index of -1). 
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Fig. 1. — The observed spectra of quasars with a central black 
hole mass of M]-,h = 10* M©. The upper curves correspond to a 
source redshift of Zs = 6 and the lower curves to a source redshift of 
2s = 11. In both cases we assume sudden reionization at Zr = 10. 
The dashed curves show the intrinsic spectral shape, taken from Elvis 
et al. (1994). We include the opacity of the neutral H and He in the 
IGM at z > Zr and the opacity due to electron scattering and the 
(extrapolated) Lya forest in the IGM at z < Zr. 

This simple model was demonstrated to accurately re- 
produce the evolution of the optical luminosity function 
in the B-band (Pei 1995) at redshifts z ^ 2.2 (Haiman 
& Loeb 1998). Because our model incorporates several 
simplifying assumptions, we regard it as the minimal toy 
model which successfully reproduces the existing data. If 
one of our input assumptions was drastically violated and 
our model had failed to fit the observed LF, then a mod- 
ification of its basic ingredients would be needed. In this 
Letter, we focus on the predictions of this minimal model 
in anticipation of the forthcoming launch of CXO; an inves- 
tigation of a broader range of plausible toy models will be 
made elsewhere. We adopt the concordance cosmology of 



Ostriker & Steinhardt (1995), namely a ACDM model with 
a tilted power spectrum (f^o, JIa, f^b, ^, fg/j-i , 7i)=(G.35, 
0.65, 0.04 0.65, 0.87, 0.96). 

Figure shows the adopted spectrum of quasars, assum- 
ing a black hole mass Mbh — IO^Mq, placed at two differ- 
ent redshifts, Zg = H and Zg = 6. In computing the inter- 
galactic absorption, we included the opacity of both hydro- 
gen and helium as well as the effect of electron scattering. 
We assumed that reionization occurred at z^ = 10 and 
that at higher redshifts the IGM was homogeneous and 
fully neutral. At lower redshifts, < z < Zr, we included 
the hydrogen opacity of the Lya forest given by Madau 
(1995), and extrapolated his fitting formulae for the evo- 
lution of the number density of absorbers beyond z = 5 
when necessary. As Figure |l| shows, the minimum black 
hole mass detectable by the ~ 2 x 10""'^^ erg s~^ cm~^ 
flux limit of CXO (see below) is Mbh ~ 10^ M© at z = 10 
and Afbh ~ 2 x 10^ Mq at z = 5. In our model, the 
corresponding halo masses are Afhaio ~ 3 x 10^^ Mq, and 
Afhaio ~ 6 X 10^° Mq, respectively. Although such massive 
halos are rare, their abundance is detectable in wide-field 
surveys. Note that an accurate determination of the spec- 
trum below ^ O.lkeV could have provided an estimate of 
the reionization redshift Zj.. Unfortunately, this spectral 
regime suffers from Galactic absorption (O'Flaherty and 
Jakobsen 1997) and is outside the 0.4-6keV detection band 
of CXO. 



3. THE X-RAY LUMINOSITY FUNCTION 

In our model, the X-ray luminosity function at a red- 
shift z [in Mpc~^ {eig/s)~^] is given by a sum over halos 
that formed just before that redshift. 



t{z) 



dt 



fxUAt) dA/bhdi 



(2) 



Mb 



-/x/q(At) 



where Lx is the observed X-ray luminosity in the instru- 
ment's detection band (0.5-3 keV for ROSAT and 0.4-6 
keV for CXO); fx is the fraction of the quasar's bolomet- 
ric luminosity emitted in this band; (PN/dMhhdt is the 
black halo formation rate, given by a convolution of the 
Press-Schechter halo mass function with equation (|l|); and 
At = t{z) — t is the time elapsed from a cosmic time t until 
a redshift z. 

Although our model was constructed so as to fit the ob- 
served optical LF, Figure || demonstrates that it is also 
in good agreement with the data on the X-ray LF. This 
implies that the choice of quasar spectrum in our model 
is reasonable. The solid curve in this figure shows the 
prediction of equation (^ at z = 3.5, near the highest 
redshift where X-ray data is available. The bottom curve 
corresponds to a cutoff in circular velocity for the host 
halos of Wcirc ^ 50 km s~^, which is introduced here in 
order not to over-predict the number of quasars in the 
HDF (Haiman, Madau & Loeb 1999). The data points 
are from recent ROSAT measurements, and the dashed 
curve in this figure shows a fitting formula from Miyaji et 
al. (1998). Note that the faintest quasar actually observed 
has Lx ^ 10"*^ erg s~^, and the fitting formula below this 
luminosity is highly uncertain. Nevertheless, it is impor- 
tant to remember that our model was calibrated so as to fit 
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the observed optical luminosity function of quasars; the ex- 
istence of a population of obscured quasars which are faint 
in the optical band but bright in X-rays could increase the 
number counts beyond our model predictions. 



shows, our quasar model predicts an unresolved flux just 
below this limit in the 0.5-3 keV range. The model also 
predicts that most ( > 90%) of this yet unresolved fraction 
arises from quasars beyond z = 5. 



z = 3.5 




Miyaji et al. 
model 
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Fig. 2. — The predicted X —ray luminosity function at z = 3.5 
in our model (solid curves). The lower curve shows the effect of a 
cutoff in circular velocity for the host halos of ii^irc ^ 50 km s~^. 
The ROSAT data points are adopted from Miyaji et al. {1998a) and 
the dashed curve shows their fitting formula (for our background 
cosmology) . 



4. THE X-RAY BACKGROUND 

Existing estimates of the X-ray background (XRB) pro- 
vide another useful check on our quasar model. The unre- 
solved background flux at a photon energy E is given by 
(Peebles 1993) 



keV ciR-^ s-^ si--^ keV-\ 



F{E) 

./n az 

(3) 

where E^ = E{1 + z)\ and j{Ez,z) is the comov- 
ing emissivity at a local photon energy Ez, in units 
of keV cm^'^ s~^ sr"-'^ keV~^, from quasars shining at a 
redshift z. This emissivity is a sum over all quasars 
whose individual observed flux at z = is below the 
ROSAT PSPC detection limit for discrete sources of 2 x 
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erg cm ^ s ^ (Hasinger & Zamorani 1997). 
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Figure || shows the predicted spectrum of the XRB in 
our model at z = (solid lines). In computing the back- 
ground spectrum, we ignored the HI absorption in the 
IGM, since it is negligible at energies above 100 eV. We 
also carried out the integral in equation (|^) only for z > 2, 
the redshift range where our model is valid (Haiman & 
Loeb 1998). The short dashed lines show the predicted 
fluxes assuming a steeper spectral slope beyond 10 keV 
{a — —0.5, or a photon index of -1.5). The long dashed 
line shows the 25% unresolved fraction of the soft XRB 
observed with ROSAT (Miyaji et al. 1998b; Fabian & 
Barcons 1992). This fraction represents the observational 
upper limit on the component of the soft XRB that could 
in principle arise from high-redshift quasars. As the figure 
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Fig. 3. — Spectrum of the soft X-ray background in our model. 
We assume that the median X-ray spectrum of each source follows 
the mean spectrum of Elvis et al. (1994) up to 10 keV, and has 
a spectral slope of 0.5 (solid lines) or -0.5 (short-dashed lines) at 
higher photon energies. The lower curves show the spectra resulting 
from quasars with redshifts between 2 < z < 5, and the upper curves 
include contributions from all redshifts z > 2. In both cases, we in- 
clude only quasars whose individual fluxes are below the ROSAT 
PSPC limit of 2 X lO"'^^ erg cm~-^ s"'^ for the detection of a dis- 
crete source. The long— dashed line shows the unresolved fraction 
(assumed to be 25%) of the soft X-ray background spectrum from 
Miyaji et al. {1998b). 



5. PREDICTED NUMBER COUNTS FOR THE CHANDRA 
X-RAY OBSERVATORY 

By summing the luminosity function over redshifts, we 
obtain the number counts of quasars per solid angle ex- 
pected to be detectable in a flux interval dFx around 
Fx = Lx /4:TTdj^{z), for all sources above a redshift z: 



dN 



dFxdn 



dz 



f dV 



V dzdQ 



<f>iLx,z)A7rdliz), (4) 



where {dV/ dzdQ) is the comoving volume element per unit 
redshift and solid angle, and dL{z) is the luminosity dis- 
tance at a redshift z. 

In Figure |, we show the predicted counts from equa- 
tion (|^) in the 0.4-6keV energy band of the CCD Imag- 
ing Spectrometer (ACIS) of CXO. Note that these curves 
are insensitive to our extrapolation of the template spec- 
trum beyond 10 keV. The figure is normalized to the 
17' X 17' field of view of the imaging chips. The solid 
curves show that of order a hundred quasars with z > 5 
are expected per field. The abundance of quasars at higher 
redshifts declines rapidly; however, a few objects per field 
are still detectable at z ~ 8. The dashed lines show the 
results for a minimum circular velocity of the host halos of 
Wcirc > 100 km s~^, and imply that the model predictions 
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for the CXO satellite are not sensitive to such a change 
in the host velocity cutoff. This is because the halos shin- 
ing at the CXO detection threshold are relatively massive, 
-^haio ~ 10^"'^ Mq, and possess a circular velocity above 
the cutoff. In principle, the number of predicted sources 
would be lower if we had assumed a steeper spectral slope. 
However, the agreement between the LF predicted by our 
model at z « 3.5 and that inferred from ROSAT obser- 
vations would be upset by such a change, and require 
a modification of the model that would in turn tend to 
counter-balance the decrease in the predicted counts. 




log {Yy^/erg s"^ cm^^) 

Fig. 4. — The total number of quasars with redshift exceeding 
z = 5, z = 7 , and 2 = 10 are shown as a function of observed X-ray 
flux in the CXO detection band. The sohd curves correspond to a 
cutoff in circular velocity for the host halos of tJcirc It. 50 km s~^, 
the dashed curves to a cutoff of ticirc > 100 km s~^. The vertical 
dashed line show the CXO sensitivity for a 5(t detection of a point 
source in an integration time of 5 X 10^ seconds. 



6. DISCUSSION 

The existence of quasars at redshifts z <; 5 has impor- 
tant consequences for the reionization history of the Uni- 
verse. Quasars can be easily distinguished from stellar sys- 
tems by their X-ray emission. We have demonstrated that 
state-of-the-art X-ray observations could provide more 
stringent constraints on quasar models than currently pro- 
vided by the Hubble Deep Field (Haiman, Madau, & Loeb 



1999). In particular, we have found that forthcoming X- 
ray observations with the CXO satellite might detect of 
order a hundred quasars per field of view in the redshift 
interval 5 z < 8. Our numerical estimates are based on 
the simplest toy model for quasar formation in a hierar- 
chical CDM cosmology, that satisfies all the current obser- 
vational constraints on the optical and X-ray luminosity 
functions of quasars. Although a more detailed analysis is 
needed in order to assess the modeling uncertainties in our 
predictions, the importance of related observational pro- 
grams with CXO is evident already from the present anal- 
ysis. Other future instruments, such as the HRC or the 
ACIS-S cameras on CXO, or the EPIC camera on XMM, 
might also be useful in searching for high-redshift quasars. 

Follow-up optical and infrared observations are needed 
in order to identify the redshifts of the faint point-like 
sources that might be detected by the CXO satellite. 
Quasars emit a broad spectrum which extends into the 
UV and includes strong emission lines, such as Lya. For 
quasars near the CXO detection threshold, the fluxes at 
~ 1/xm are expected to be relatively high, ^ 0.5-0.8 ^Jy. 
Therefore, infrared spectroscopy of X-ray selected quasars 
could prove to be a particularly useful approach for unrav- 
eling the reionization history of the intergalactic medium 
at 2 ^ 5 (for the potential lessons to be learned, see 
Miralda-Escude 1998; Haiman & Loeb 1999a; Loeb & Ry- 
bicki 1999; and Rybicki & Loeb 1999). 

At present, the best constraints on hierarchical models 
of the formation and evolution of quasars originate from 
the Hubble Deep Field. However, HST observations are 
only sensitive to a limiting magnitude oi V and 
cannot probe the earliest quasars. The Next Generation 
Space Telescope (NGST), scheduled for launch in 2008, 
will achieve nJy sensitivity at wavelengths 1 — 3 ^m, and 
could directly probe the earliest quasars (Haiman & Loeb 
1998). The combination of X-ray data from the CXO 
satellite and infrared spectroscopy from NGST could po- 
tentially resolve one of the most important open questions 
about the thermal history of the Universe, namely whether 
the intergalactic medium was reionized by stars or by ac- 
creting black holes. 
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and the NASA grant NAG 5-7092 at Fermilab. AL was 
supported in part by NASA grants NAG 5-7039 and NAG 
5-7768. 
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